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Abstract 

The program SOFTSUSY can calculate tree-level neutrino masses in the R -parity violating minimal supersymmetric 
standard model (MSSM) with real couplings. At tree-level, only one neutrino acquires a mass, in contradiction 
with neutrino oscillation data. Here, we describe an extension to the SOFTSUSY program which includes one-loop 
/^-parity violating effects' contributions to neutrino masses and mixing. Including the one-loop effects refines the 
radiative electroweak symmetry breaking calculation, and may result in up to three massive, mixed neutrinos. This 
paper serves as a manual to the neutrino mass prediction mode of the program, detailing the approximations and 
conventions used. 
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1. Program Summary 

Program title: SOFTSUSY 

Program obtainable from: http : //projects . hepf orge . org/ sof tsusy/ 



Distribution format: tar.gz 
Programming language: C++, f ortran 
Computer: Personal computer 
Operating system: Tested on Linux 4.x 
Word size: 32 bits 
External routines: None 

Typical running time: A second per parameter point. 

Nature of problem: Calculation of neutrino masses and the neutrino mixing matrix at one-loop level in the /^-parity 
violating minimal supersymmetric standard model. The solution to the renormalisation group equations must be con- 
sistent with a high or weak-scale boundary condition on supersymmetry breaking parameters and /^-parity violating 
parameters, as well as a weak-scale boundary condition on gauge couplings, Yukawa couplings and the Higgs poten- 
tial parameters. 

Solution method: Nested iterative algorithm. 

Restrictions: SOFTSUSY will provide a solution only in the perturbative regime and it assumes that all couplings of 
the MSSM are real (i.e. CP-conserving). 
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2. Introduction 



Supersymmetric (SUSY) models provide an attractive weak-scale extension to the Standard Model (SM). The 
R-parity conserving (R p ) minimal supersymmetric extension of the Standard Model (R p MSSM) is often used as a 
reference model for phenomenological studies. There exist several publicly available spectrum generators for the R p 
MSSM: ISASUGRA SOFTSUSY [2], SUSPECT Jl and SPHENC0 Jj. Spectrum information is typically transferred 
to decay packages and event generators via a file in the SUSY Les Houches Accord format JHH]. 

The most general renormalisable superpotential of the MSSM contains /^-Parity violating (R p ) couplings which 
violate baryon and lepton number 10]. A symmetry can be imposed upon the model in order to maintain stability 
of the proton, for example baryon triality Jll] or proton hexality 131 . It has been shown that Ijtp models may have 
interesting features, such as the generation of non-zero neutrino masses and mixing without the addition of right- 
handed neutrino fields ifioll . In fact, neutrino oscillation data indicates that at least two neutrinos must be massive, so 
a realistic extension to the SM should include mechanisms to generate these masses. Here, we describe an extension 
to SOFTSUSY which calculates neutrino masses and mixing to one-loop order in the presence of Ijt p couplings. The 
latest version of SOFTSUSY including Ijt p effects can be downloaded from the address 



http: //projects . hepf orge . org/sof tsusy/ 



Installation instructions and more detailed technical documentation of the code may also be found there. 

The R p and R p aspects of the SOFTSUSY calculation leading to self-consistent spectra are already explained in 
detail in Refs. ]2|] and 111 111 respectively, with the technical differences between these two calculations detailed in 



the latter reference. They shall not be repeated here. Instead we shall concentrate on the calculation of neutrino 
masses and mixing, and differences in the ft p calculations between the new and the previous release, which are called 
the neutrino mode and the Ijt p mode respectively. Two new features of S0FTSUSY3 . 2 are the tji p one-loop tadpole 
corrections to the two Higgs vacuum expectation values (VEVs), and the complete one-loop tadpole corrections to 
the sneutrino VEVs. These improvements allow the computation of the neutrino spectrum and neutrino mixings, with 
minimal additional computational cost. They also slightly change some predictions in the non-neutrino sector. 

We proceed with a definition of the SOFTSUSY convention for the Ijt p parameters and mixings relevant for the 
neutrino mass calculation in section[3] Next, in section|4] we discuss the calculation of the neutrino masses, detailing 
the approximations made. Installation instructions can be found on the SOFTSUSY web-site, but instructions to run the 
program can be found in|A] The output from a SOFTSUSY sample run in the neutrino mode is displayed and discussed 
inlBl whereas a sample main program is shown and explained in|C] Some more technical information on the structure 
of the program can be found in|D] It is expected that the information in [D] will only be of use to users who wish to 
'hack' SOFTSUSY in some fashion. 



3. MSSM ft p Parameters 

In this section, we introduce the Ijt p MSSM parameters in the SOFTSUSY conventions relevant for the neutrino mass 
and 1-loop tadpole calculations. A detailed description of the complete set of Iji p MSSM parameters is presented in 
Ref. Illll . The latter follows Ref. IU2I1 and so the notation and conventions employed are similar. 

3.1. Iji p supersymmetric and SUSY breaking parameters 

The chiral superfield particle content of the MSSM has the following S U(3) c xS U(2)lX U(1)y quantum numbers: 

L: (1,2,-1), E - (U.D. G:(3,2,i), U:(3,l,-h, 

6 3 

D: (3,1, |), Hi: (1,2, -i), H 2 : {\,2,\). (D 



SPHENO includes a small subset of the^?,, interactions. 
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L, Q, Hi, and Ho are the left handed doublet lepton and quark superfields and the two Higgs doublets. E, U, and D 
are the lepton, up-type quark and down-type quark right-handed superfield singlets, respectively. The ft p part of the 
renormalisable MSSM superpotential that violates lepton number, written in the interaction eigenbasis, is 



-A ijk L a ,L)E k + A' ijk LfQfD kx - k^H* 



(2) 



Here, k, [Atj k and A' rk ] are the bi-linear [trilinear] couplings. We denote an S U(3) colour index of the fundamental 
representation by e {1,2,3}. The S U{2)l fundamental representation indices are denoted by \a,b,c] e {1,2} 

and the generation indices by [i, j, k] € {1,2,3}. e a i, = ^ b is a totally antisymmetric tensor, with = 1. Currently, 
only real couplings in the superpotential and Lagrangian are included. 

The corresponding soft JjL p breaking interaction potential, together with the bi-linear mixing term between the 
scalar component of the lepton doublets and the H\ superfields, is involved in tadpole and 1-loop corrections to the 
neutrino-neutralino masses. It is given by 



1 



^h i j k L i i 'L b j ~e k + h' ijk L^Qfd kx - DiLfH* + H.c 



bxj 



\ml „ V. ff? + H.c.l 

L L,H\ ia 1 J 



(3) 



where fields with a tilde are the scalar components of the superfield with the same capital letter. The scalar components 
of the superfields Hi and H2 and the superfields themselves have the same notation. The electric charges of d and e 
are |, and 1, respectively, ftyj and h'^ k are trilinear soft SUSY breaking parameters that correspond to Ajj k and A'.-, 
and Dj and m? are bi-linear SUSY breaking parameters. "H.c." denotes the Hermitian conjugate of the preceding 
terms. 



3.2. Neutral and charged fermion masses 

In the presence of R p interactions that violate lepton number, the neutrinos and neutralinos mix, and the charginos 
and leptons mix with each other. Sneutrino-anti-sneutrino mixing is also present in principle; in practice this has been 
shown to have negligible phenomenological consequences once experimental bounds have been applied II 1311 . and is 
neglected in our calculation. 



The (7x7) neutrino-neutralino Lagrangian mass term, containing three families of neutrinos is given in II 1 211 and 
reads 



£ = -- (v h -m, - w (3) , E, h° 2 ) m n 



(4) 



where at tree level, 
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(5) 



In Eq. (0), the flavour basis neutral fermions are the 3 neutrinos (v,), the gauginos {HfW) and the higgsinos {h Q v h®). 
In the mass matrix in Eq. (0, k-, [yu] are the supersymmetric Ijtp [R p ] bi-linear mixing parameters, v, [(Hp and (H®)] 
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are the sneutrino [H 2 and H®] vacuum expectation values ( VEVs), Mi , M 2 are the gaugino masses of hypercharge and 
weak isospin respectively, and g', g 2 are the corresponding gauge couplings. The mass eigenstates are obtained upon 
diagonalisation of Ai: ^=1,2,3, x\ 234 v * a a 7 by 7 orthogonal matrix O: 



At 



diag 



T M N o. 



(6) 



A simple multiplication of rows of O by factors of / can absorb any minus signs in M^" 8 . Mass eigenstates are ordered 
in increasing (absolute) mass eigenvalues. 

To facilitate comparisons with neutrino oscillation data, it is useful to express AIn as 



Ms = 



m v m 
m T Mu 



(7) 



where m v is a 3 x 3 mass matrix, Ai^o is a 4 x 4 mass matrix, and m is a 3 x 4 matrix that mixes the neutrinos and the 
neutralinos. We define an effective 3x3 neutrino mass matrix A4 V via the see-saw relation 



M v = m v - m Mj m T . 



It can be diagonalised by a (9(3) matrix O v , 



M d J ag = OlM V Oy. 



(8) 
(9) 



For realistic mass spectra, O v is practically the same as the 3x3 neutrino part of O (we do not assume this, however). 
At tree level, the effective neutrino mass matrix is 112 ] 



AI" 



fi(M lg 2 2 + M 2 g' 2 ) 



2<// ( 1 ) >(//?>(M I ^2 + M 2 g' 2 ) - 2nM x M 2 



A1A1 AiA 2 A1A3 
A 2 Ai A 2 A 2 A 2 A 3 
A3A1 A3A2 A3A3 J 



where 



A, = v, --(Hi)-, 



i= 1,2,3. 



(10) 



(11) 



A4y ee is rank 1 and so contains two zero eigenvalues. Because of the presence of two massless neutrinos in the 
spectrum, the tree level values are not realistic. The numerical values of O v are also not meaningful in the tree-level 
approximation. The 1-loop corrections to A4 V can lead to two or three non-zero neutrino masses, if more than one 
lepton flavour is violated. 

By default, the neutrino masses are normal-ordered, i.e. \m Yl \ < \m Vl \ < |m V3 |. A spectrum with inverted ordering 
can be obtained by 
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together with the corresponding swaps (of the column vectors) in the mixing matrix O v . The mass ordering then 
becomes |ra v <| < |m v 'J < |m v ;|. 

The presence of Ijt p interactions also mix charged-leptons with the charginos. The Lagrangian contains the (5 X 5) 
chargino-lepton mass matrix 



£ = -(-iW-,h- l ,e- L .)M c 



e R k ) 



+ H.c. 



(13) 



The mass eigenstates I = (e,//, t),x? 2 are given upon the diagonalisation of the matrix Ale, where at tree level, 



M c 



I M 2 f 2 {H° 2 ) 0j 

f 2 Vi Kt X ((Y E ) U (H1) + X UjVk ) J 



(14) 



Here, Ye is the lepton Yukawa matrix from the R p superpotential in Ref. We define the diagonalised mass matrix 

Mf g = UM C V T , (15) 

U and V being orthogonal 5x5 matrices. The one-loop fji p corrections to Mc are related to neutrino masses. For# ; , 
effects giving small enough neutrino masses to pass empirical bounds, the one-loop Ijt p corrections to Aic are expected 
to be negligible, unless one assumes that the Ijt p effects are large, but have a high degree of cancellation in the neutrino 
masses. At present, the one-loop Ijt p corrections to Aic are not implemented. The absolute mass eigenvalues are in 
increasing order along the diagonal of M Q mg '. An effective charged lepton mixing matrix, Uu can be obtained from 
the 3x3, top right part of U. The PMNS matrix t/pMNS is then defined as 

[/pmns = u;o v , (16) 

with O v including the (column) swapping when an inverted mass ordering is desired. 



3.3. Neutral slepton masses 

Lepton number violating interactions mix the sleptons and higgs'. In the limit of CP conservation, as assumed 
here, one obtains two mass squared matrices: one for the neutral CP even real scalars At^ , and one for neutral CP 

odd real scalars Ai 2 . The Lagrangian containing these terms looks like 



( Vk± J 



(17) 



where the fields with subscript +(-) are the real (imaginary) part of the complex fields in obvious notation. At tree 
level, 



{H° } 2 -<#°><#°> -<#>* ) 
-(H° 2 ) Vj (H°) Vj v jVk 



"2 



+m 



-Dj fiKj + m z L Hj KjK k + (mp jk ) 



( s 
1 

Sjk , 



(18) 



where m\ is the soft breaking term corresponding to the bilinear Higgs-mixing parameter //. 

Currently, yVP is calculated at tree-level and is only used to calculate the one-loop corrections to neutrino masses 
and mixings. The mass eigenstates are ordered in mass (m o) i=1 5 



diagC/n^ = J ± h?^0 ± 



(19) 



with mixing matrix ± . 



4. Calculation Algorithm 

The calculation to calculate the low energy ft p mass spectrum proceeds mostly in the same way as the original 
Iji p mode, which is detailed in Ref. Hill , but now also includes the full 1-loop Ijt p tadpoles and additional methods 
to calculate the 1-loop 7x7 neutralino-neutrino mass matrix and to obtain the neutrino pole masses and the PMNS 
matrix via the seesaw-mechanism. The calculation algorithm is displayed in Fig.Q] where the changes with respect to 
the Ijt p mode are displayed in bold print. For a given a set of inputs forR p couplings, assumptions about the measurable 
quark mixing make a physical difference to the results. It is well known that in the R p MSSM, a t/(3) 5 global family 
symmetry, where Yukawa couplings also transform as spurions, renders physical results invariant to where the CKM 
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1. SUSY radiative corrections to gi(Mz) 



2. Ye, Yq include va contributions. Iterative solution of Ye- 



convergence 



3. Run to M s . 



4. EWSB, iterative solution of p and sneutrino VEVs 
including full R p and R p 1 loop contributions to the tadpoles jj— - 



5. Run to Mx- Apply soft breaking and# ; , SUSY boundary conditions. 



6. Run to M z . 



7. Calculate Higgs, sparticle and neutrino pole masses. Run to Mz- 



Figure 1: Iterative algorithm used to calculate the Ijtp MSSM spectrum in the neutrino mode. Differences to the original ft„ mode are highlighted 
in bold print. The initial step is the uppermost one. Mg is the scale at which the EWSB conditions are imposed, as discussed in the text. Mx is the 
scale at which the high energy SUSY breaking boundary conditions are imposed. 



quark mixing lies. For example, the physics is identical for the cases where CKM mixing lies solely in left handed 
up quarks and for the case where it lies solely in left handed down quarks, since the family symmetry relates the 



two cases 111 411 . However, in the ft p MSSM, one would also have to transform the Ijtp Yukawa couplings A!^ k and A"j k 
to keep the theory invariant under the transformation. Since this is inconvenient, we wish to leave the Ijt p couplings 
untransformed, and instead specify (in the basis in which the Ijt p couplings are given), where the quark mixing lies. 
There is thus a physical difference depending on whether one assumes that CKM mixing lies e.g. solely in the down- 
type quarks or solely in the up-type quarks. We describe the parameter in the program which specifies this mixing 
inlDl It can have a significant effect on the neutrino mass predictions. 

4.1. Ijt p tadpole contributions 

The 1-loop ]jt p contributions to the tadpoles 1^, as well as R p and Jjt p contributions to are implemented in 
the neutrino mode of SOFTSUSY . These calculations include all three SM fermion families in the loop. In the presence 
of Ijt p interactions, the charged sleptons and down-type Higgs mix, and the neutral scalars and Higgs mix with each 
other, so out computation treats these contributions on an equal footing. The sneutrino VEVs are now also calculated 
at 1-loop level at the renormalisation scale Ms = a/»z,~ (Ms ) »z,~ (Ms ). However,there are no renormalisation group 
equations (RGEs) implemented for the sneutrino VEVs, so they remain constant at any renormalisation scale. 

4.2. Neutrino-neutralino mass corrections 

The physical neutrino-neutralino masses are calculated to full one-loop order, including both R p and Ijt p contribu- 
tions, following Refs. Il5lll6ll closely. The running parameters are evaluated at the renormalisation scale p = Msusy- 
However, note that the physical neutralino masses are displayed in the R p limit in the SOFTSUSY output without in- 
cluding the ft p contributions. However, in the absence of large cancellations, Ijtp corrections to neutralino masses will 
be of order the neutrino mass, i.e. negligible. To calculate the physical neutrino masses, the external momentum for 
the physical neutrino-neutralino mass matrix, M^ ys , is set at p = 0. The matrix M'^ >s is given by 

A< ,v ' 5 = M N + i(A N + A£), (20) 
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where 

A N = E D -M N I L . (21) 

In the above expression, "Zq and are mass corrections and wavefunction renormalisation respectively. The effective 
physical neutrino mass matrix At v is obtained as defined in Eq. [8] from which the physical neutrino masses and 
mixings are extracted. 

To avoid numerical instability due to large cancellations between the CP-even (CPE) and CP-odd (CPO) Y.d 
contributions to m v in Eq. (0, the combined contribution is obtained using an analytic expansion in the 2 x 3 ft p 
matrices that mix the Higgs with the sneutrinos in the 5x5 CPE and CPO neutral scalar mass matrices. Both R p 
and ft p corrections to the physical neutrino masses are computed simultaneously, treating the charged sleptons and 
down-type Higgs as well as the neutral scalars and Higgs particles on an equal footing, as in the tadpole calculation. 

When calculating realistic inverted or quasi-degenerate neutrino mass spectra, the mass eigenvalues are much 
larger than their differences, so some fine tuning is expected. For example, in Ref. Il7ll . the level of fine-tuning is 
estimated to be of O(10~ 3 ) and <9(1(T 4 ) for the inverted and quasi-degenerate cases respectively. In order to obtain 
numerically stable results for neutrino masses and mixings for degenerate or inverted hierarchies, we recommend 
setting the TOLERANCE parameter to 1(T 5 and 1(T 6 , respectively. 
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A. Running SOFTSUSY 

SOFTSUSY produces an executable called sof tpoint .x. Instructions on how to run sof tpoint . xin the presence 
of ftp couplings are detailed in Appendices A and B of Ref. Hill . Starting from version SOFTSUSY-3 . 2, output for 
both input options is compliant to the SLHA2 conventions, which are fully explained in Ref. [6]. A sample SLHA2 
file, rpvHouchesInput, may be used as input, and SOFTSUSY run, by the command 

. /sof tpoint . x leshouches < rpvHouchesInput 

When any non-zero ft p couplings are set in either the command line or the ft p Les Houches input file, the program 
automatically calls the neutrino mode instead of the original ft p mode. There is a new optional SLHA2 input parameter 
in the block SOFTSUSY associated with the neutrino mode, which can be set in the following fashion: 

Block SOFTSUSY 

9 0.000000000e+00 # output uses normal hierarchy (=0.0) or inverted (=1.0) 

We also remind the reader that the quark mixing assumption can have a physical effect given an assumed set of ft p 
couplings, and that the CKM mixing may be set solely in the up or down quark sector via 

Block SOFTSUSY 

2 1.000000000e+00 # quark mixing: none (=0), up (=1) or down (=2) 

B. Sample Output 

The following is an example for one non-zero ft p coupling A'^ 2 {Mcut) = 0.01 at Mgut in an ft p mSUGRA 
SPS la-like JH] scenario: 

. /sof tpoint .x sugra 100 250 -100 10 unified 1 lambdaP 3 2 2 0.01 

The SLHA2 compliant output contains /^,-conserving and ft p MSSM information. Results related to the neutrino 
sector are displayed in the following blocks: 
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Block MASS # Mass spectrum 

#PDG code mass particle 

12 4. 128113522451627e-50 # Mnul(pole) normal hierarchy output 

14 -2.818848508640361e-42 # Mnu2(pole) normal hierarchy output 

16 -5.258228488951831e-06 # Mnu3(pole) normal hierarchy output 

Block RVSNVEV Q= 4 . 670086691357357e+02 # sneutrino VEVs D 

1 -2. 115552831665844e-23 # SneutrinoVev_{l> 

2 -3. 117652325269568e-18 # SneutrinoVev_{2> 

3 -7.695902798823816e-02 # SneutrinoVev_{3> 



Block UPMNS Q= 4 . 670086691357357e+02 # neutrino mixing matrix: 

# UPMNS_{11> matrix element 

# UPMNS_{12> matrix element 

# UPMNS_{13> matrix element 

# UPMNS_{21> matrix element 

# UPMNS_{22> matrix element 

# UPMNS_{23> matrix element 

# UPMNS_{31> matrix element 

# UPMNS_{32> matrix element 

# UPMNS_{33> matrix element 



9 . 999999999518719e-01 
9.811018586559509e-06 
2 . 652785268618422e-20 
-9.811018586559507e-06 
9 . 999999999518718e-01 
6 . 939913220612296e-17 
-2.454643369143216e-20 
-2 . 018886379347570e-16 
1 . 000000000000000e+00 



The block MASS includes the three pole neutrino masses calculated in units of GeV, given the ]jt p model parameters. 
Note that the corresponding particle data group (PDG) codes lfl9ll are for the neutrino flavour eigenstates. However, 
since we are dealing with massive neutrinos, we slightly adapt this definition such that the mass ordering of the neu- 
trinos is specified to be either normal or inverted hierarchy according to the value of the invertedOutput parameter. 
In block RVSNVEV, the 1-loop sneutrino VEVs are displayed at the energy scale Q at which the EWSB conditions 
are imposed (Ms). Finally, the block UPMNS displays the matrix elements of the matrix I/pmns m Eq.Q~6] The neu- 
trino spectrum predicted by this model point is in contradiction with data, but in the next section, we calculate with a 
parameter point that leads to neutrino mixings and masses that are agreement with 201 1 oscillation data in a sample 
program. 



C. Sample Program 

In this section we present a sample main program, that illustrates a calculation of neutrino masses and mixings 
given a set of Ijt p parameters at M<jut ■ This main program is included as the rpvneutmain . cpp file with the standard 
SOFTSUSY distribution and performs the calculation assuming mSUGRA parameters m = 100GeV,Mi/ 2 = 500 GeV, 
A = 932.9 GeV, tan/? = 20 and fi > 0, with non-zero ljt p parameters A' m (M GUT ) = 0.039, A' 2n (M CUT ) = -0.016, 
A' 3n (M CUT ) = 0.018, A\ 33 (M CUT ) = 3.0 x lO" 5 , A' 233 (M CUT ) = 3.0 x 10" 5 , A' 333 (M CUT ) = -3.5 x lO" 5 . 

The sample program has the following form: 

#include <rpvmain.h> 

int main() { 

/// Sets up exception handling 
signal (SIGFPE, FPE_ExceptionHandler) ; 

/// MIXING=0: diagonal Yukawa matrices, MIXING=1: CKM-mixing in up-sector, 
/// MIXING=2: CKM-mixing in down-sector 
MIXING = 1; 

/// Apply SUSY breaking conditions at GUT scale, where g_l=g_2 
bool gaugeUnif ication = true; 
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/// Sets format of output: 3 decimal places 
outputCharacteristics (3) ; 



/// "try" catches errors in main program and prints them out 
try { 

QedQcd oneset; ///< See "lowe.h" for default parameter definitions 

oneset .toMz() ; ///< Runs SM fermion masses to MZ 

/// Guess at GUT scale 
double mxGuess = 2.el6; 

/// Close to scenario IH S2 from arXiv 1106.4338 
int sgnMu = 1 ; 

double tanb =20., aO = 912.3, ml2 = 500., mO = 100.; 

/// Define RpvNeutrino object 
RpvNeutrino kw; 

/// Set the GUT scale RPV SUSY couplings 
kw.setLamPrimed, 1, 1, 0.040); 
kw.setLamPrime(2, 1, 1, -0.018); 
kw.setLamPrime(3, 1, 1, 0.019); 
kw.setLamPrimed, 3, 3, 3.3e-5); 
kw. setLamPrime (2 , 3, 3, 3.2e-5); 
kw. setLamPrime (3, 3, 3, -3.5e-5); 

/// Store inputs into one vector 

DoubleVector pars(3); pars(l) = mO; pars(2) = ml2; pars(3) = aO; 

/// Outputs the RPV couplings required into the vector pars used by lowOrg 
kw.rpvDisplay(pars) ; 

/// Makes sure the neutrino mass ordering will be as expected in inverted 
/// hierarchy output. If required, must be set before lowOrg is called 
kw. setlnvertedOutput ; 

/// For inverted output, because of large cancellations, one requires the 
/// SUSY spectrum calculation to be performed to a high fractional 
/// accuracy 
TOLERANCE = 1.0e-6; 

/// Main driver routine: do the calculation 

double mgut = kw. lowOrg (rpvSugraBcs , mxGuess, pars, sgnMu, tanb, oneset, gaugeUnif ication) ; 
/// Output the results in SLHA2 format 

double qMax =0.; char * modelldent = (char *)"sugra"; 
int numPoints = 1; bool altEwsb = false; 

kw. lesHouchesAccordOutput (cout , modelldent, pars, sgnMu, tanb, qMax, numPoints, mgut, altEwsb); 

} 

catch(const string & a) -[ 
cout « a; exit(-l); 

} 

catch(const char *a) { 

printf 07.S" , a); exit(-l); 

} 

} 

The structure of the main file above is as follows. After including a header file, an exception handler is set up, fol- 
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lowed by specifying a (sub-) set of global variables, all of which are described in the R p manual J2l- These include the 
MIXING switch, which determines how any quark mixing is implemented, the gaugeUnif i cat ion switch, which de- 
termines mGutGuess as the scale of Mqut of electroweak gauge unification. The switch outputCharacteristics, 
which specifies the output accuracy, is then set. 

The running masses of the SM fermions and the QED and QCD gauge couplings are determined at Mz from data 
with the method toMz. In order for SOFTSUSY to determine Mcut, it requires an initial guess, which must be supplied 
as the initial value of the variable mxGuess (in GeV), and is later over-written by the program with a more accurate 
calculated valu^l 

Next, the mSUGRA parameters sgn// =sgnMu, tan/3 =tanb, Ao/GeV=aO, Mi/2/GeV=ml2 and mo/GeV=mO are 
defined. This is then followed by the instantiation of the RpvNeutrino object kw and the assignment of the Ijtp 
parameters kw.setLamPrime(i, j ,k,val) for X { . k =val. In the iterative SOFTSUSY algorithm the parameters in 
the RpvNeutrino object change due to the RGE running. The pars vector is needed to keep track of the bound- 
ary conditions set at Mqut ■ These boundary conditions are re-set in every iteration at M<jut from the unchanged 
DoubleVector pars parameters. 

We do not fill the other 102 ft p entries of pars explicitly. This would be tedious and an additional source of 
potential bugs. Instead, we fill the RpvSof tsusy object itself using the setLamPrime method in this example. We use 
the rpvDisplay method: this fills the pars vector automatically with what was set already inside the RpvNeutrino 
object, while leaving the first nine entries in the vector unchanged. The rpvDisplay method automatically changes 
the length of pars appropriately. 

The neutrino masses may be presented in normal (default) and inverted mass orderings. This is controlled by the 
method setlnvertedOutput which must be invoked before the actual SOFTSUSY main driving method lowOrg is 
called. When lowOrg is called, the first argument specifies the type of boundary condition (currently rpvSugraBcs), 
which assumes that pars has already been prepared by using the rpvDisplay object. This is followed by a print to 
standard output of the model parameters at Mz and the physical parameters SLHA2 compliant format. Finally, the 
catch commands print any errors produced by the code. 

The output is in standard SLHA2 format, including neutrino masses in GeV units, neutrino mass ordering (inverted 
hierarchy in this example) and the PMNS mixing matrix: 

[ ... ] 

Block MASS # Mass spectrum 

# PDG code mass particle 
[ ... ] 

12 4.65793809e-ll # Mnul inverted hierarchy output 
14 4.70348794e-ll # Mnu2 inverted hierarchy output 
16 1 .20158889e-14 # Mnu3 inverted hierarchy output 

[ ... ] 



.ock UPMNS 1 


5= 9.11876000e+01 


# neutrino 


mixing 


matrix (inverted 


1 


1 


8 


.61178721e 


-01 


# 


UPMNS. 


.{11} 


matrix 


element 


1 


2 


5 


.08244716e' 


-01 


# 


UPMNS. 


.{12} 


matrix 


element 


1 


3 


7 


.64979577e- 


-03 


# 


UPMNS. 


.{13} 


matrix 


element 


2 


1 


-3 


.29127098e 


-01 


# 


UPMNS. 


.{21} 


matrix 


element 


2 


2 


5 


.69021210e 


-01 


# 


UPMNS. 


.{22} 


matrix 


element 


2 


3 


-7 


.53584909e' 


-01 


# 


UPMNS. 


.{23} 


matrix 


element 


3 


1 


-3 


.87358444e' 


-01 


# 


UPMNS. 


.{31} 


matrix 


element 


3 


2 


6 


.46453533e' 


-01 


# 


UPMNS. 


.{32} 


matrix 


element 


3 


3 


6 


.57306066e' 


-01 


# 


UPMNS. 


.{33} 


matrix 


element 
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RGE 




QedQcd 


MssmSusy 






FlavourMssmSoftsusy 


SoftParsMssm 












RpvSof tsusy 




MssmSoftsusy 



RpvSof tPars 



RpvSusyPars 



RpvNeutrino 



Figure 2: Heuristic high-level object structure of SOFTSUSY-3.2. Inheritance is displayed by the direction of the arrows. 

D. The RpvNeutrino class 

We now go on to sketch the class RpvNeutrino. This class publicly inherits from the class RpvSof tsusy, 
cf. Fig. The RpvSof t susy methods physical, calculateSneutrinoVevs and doCalcTadpoleloneLoop, 
doCalcTadpole2oneLoop are overloaded in RpvNeutrino in order to include 1-loop neutrino masses, sneutrino 
VEVs and tadpoles. The data and methods in the RpvNeutrino class deemed of possible importance for prospective 
users are presented in TableQ] Note that different quark mixing assumptions can significantly affect the predictions of 
the neutrino masses. The MIXING parameter is implemented in the neutrino mode in the same fashion as in Ref. J2l, 
and can take values , 1 , 2 for no mixing, up quark mixing and down quark mixing, respectively. 
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data variable 



methods 



DoubleVector 
"Vj=iA3 /GreV 


physNuMasses 


physical neutrino masses 


display PhysNuMasses 
setPhysNuMasses 


DoubleMatrix 
t/pMNS (3 X 3) 


uPmns 


PMNS mixing matrix 


d i s p 1 ay Upmn s 
setUpmns 


DoubleMatrix 
0(7x1) 


physNeutMix 


physical neutral fermion mixing 


displayPhysNeutMix 
setPhysNeutMix 


bool invertedOutput 


neutrino mass ordering 


display InvertedOutput 
set InvertedOutput 
setNormalOutput 


double thetal2, theta23 
double thetal3, deltaCkm 
On, 023 , 0\3, <5/radians 


PDG parameterisation 
of CKM angles 


display ThetaCkml2 , display ThetaCkm23 
display ThetaCkml3 , displayDeltaCkm 
setThetaCkml2, setThetaCkm23 
setThetaCkml3 , setDeltaCkm 


DoubleVector 

{m lf p) i= i...5 


CPEmasses 


tree-level CP even neutral 
scalar masses 


displayCPEMasses 
setCPEMasses 


DoubleVector 

(»y>=i...5 


CPOmasses 


tree-level CP odd neutral 
scalar masses 


displayCPOMasses 
setCPOMasses 


DoubleMatrix 

M 2 B (5 x 5) 

v+ 


CPEscalars 


tree-level CP even neutral 
scalar mass 2 matrix 


displayCPEscalars 
setCPEscalars 


DoubleMatrix 
M 2 (5x5) 


CPOscalars 


tree-level CP odd neutral 
scalar mass 2 matrix 


displayCPOscalars 
setCPOscalars 


DoubleMatrix 
+ (5 x 5) 


CPEscalarMixing 


tree-level CP even neutral 
scalar mixing matrix 


displayCPEscalarMixing 
setCPEscalarMixing 


DoubleMatrix 
(5 x 5) 


CPOscalarMixing 


tree-level CP odd neutral 
scalar mixing matrix 


displayCPOscalarMixing 
setCPOscalarMixing 



Table 1: RpvNeutrino class. We display the important data contained within the object, along with ways of accessing and setting their values. 
For the definitions of the methods, see the file rpvneut .h in the SOFTSUSY distribution. 
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